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Abstract. - Recent thermal conductivity measurements in the heavy-fermion compound Celrlns 
indicate that its superconducting order parameter is very different from CeCoIns. Here we show 
that these experiments are consistent with chiral d-wave symmetry, i.e. A(k) ~ e*"* cos(cfc z ). 



The discovery of antiparamagnon mediated superconductivity in the 115 heavy-fermion 
compounds CeTIn 5 , where T represents Co, Ir, Rh, or a mixture of these, has recently 
opened up a new avenue to unconventional nodal superconductivity [1]. These strongly 
interacting materials are characterized by a plethora of competing ground states in addition 
to superconductivity, including conventional and unconventional spin density wave (SDW) 
phases [2]. Among the 115 compounds, the currently most well studied is CeCoIn 5 for which 
a d-wave superconducting order parameter A(fc) <~ cos(2</>) = k x — ky has been identified 
[3-6]. Indeed, there are many parallels between CeCoIns and the high-T c cuprates, including 
(a) a layered quasi-two-dimensional Fermi surfaces [7], (b) d-wave superconductivity, and 
(c) d-wave spin density wave order in the pseudogap phase [8-10]. 

Recent thermal conductivity measurements [11, 12] indicating an order parameter sym- 
metry in Celrlns very different from the one in CeCoIns came as a big surprise. An initial 
analysis of this data suggested a hybrid E g gap, A(fc) <~ ^2,±i(^) <fi), based on the assump- 
tion that the Fermi surface is three-dimensional. However, the Fermi surface of Celrln 5 is 
in fact quasi- two-dimensional, as known from band structure analysis [7,12]. Therefore, one 
needs to consider instead superconductivity in layered structures, similar as discussed in 
Refs. [13, 14]. In this case, only / = e ±l *sin(x) (chiral d-wave) with \ — °kz or / ~ sin(x) 
(non-chiral p-wave) are consistent with the observed thermal conductivity data [11]. The 
magnitudes of the d-wave and chiral d-wave/non-chiral p-wave order parameters |A(fc)| are 
shown in Fig. 1. 

In the following, we present a theoretical analysis based on a generalized BCS model 
that properly accounts for a quasi-two-dimensional Fermi surface and a chiral d-wave super- 
conducting order parameter. The thermal conductivity is computed following the receipe 
given in Refs. [14,15]. Here, we assume for simplicity that the quasiparticle scattering is 
due to impurities. Furthermore, we consider the physically relevant limit T/A <C 1, where 
r is the quasiparticle scattering rate in the normal state and A(= 0.856K) is the maxi- 
mum value of the energy gap at T = OK. This A is the weak-coupling value for nodal 
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Fig. 1: Magnitude of order parameters |A(fe)| for d-wave (left) and chiral d-wave/non-chiral p-wave 
(right) superconductors. 



superconductors [13,16]. 

Let us begin by considering the zero-temperature limit. For quasi- two-dimensional 
structures, the thermal conductivity strongly depends on the direction within the mate- 
rial. Therefore we will discuss the cases q\\a (in-plane) and q\\c (out-of-plane) separately. 
Here q denotes the heat current. For q\\a, one obtains 
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and similarly for <f||c 
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where r a and T c denote the in-plane and out-of-plane scattering rates respectively. Eq.l 
describes the universal heat conduction as discovered by P. Lee [17,18], whereas Eq.2 is 
very different. The strength of the impurity scattering can be extracted directly from the 

r 

experimental data show in Fig. 2 of Ref. [11], from which we can deduce that — = 0.19635. 
Furthermore, from the observed anisotropy of the thermal conductivity, we can infer the 



ratio of the Fermi velocities along the c-axis and the a-b plane, i.e. 



0.66, which is very 



similar to — = 0.5 extracted for CeCoIns 
the regime T>r, 



Then, for T ^ OK but — < 1, we obtain in 
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This is consistent with the experimental observation of a dominant in-plane heat conduc- 
tivity propertional to the temperature, and a subdominant out-of-plane conductivity. 

In order to connect these finite-temperature results with the above equations for T = 0, 
we use an interpolation formula which applies in the regime for T/A(T) <ti 1. The resulting 
low-temperature thermal conductivities are then given by 
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Fig. 2: Thermal conductivity in the q\\a (in-plane) and <f||c (out-of-plane) directions. T c — OAK. 
The symbols represent experimental data from [11], and the solid lines are low-temperature fits 
using Eqs. 5 and 6. 
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respectively. 

In Fig. 2, we compare these dependencies with the experimental data reported in 

r c 

Ref. [11]. A fit of the low-temperature regimes yields good agreement with — - = 0.5592. 
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Evidently, the quasi-particle scattering rate is somewhat anisotropic in the present system. 
Here, the temperature dependence of the gap function A(T), is approximated by 



A(T) = 2.14T C 
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with T c = OAK , which is known to be a very good approximation for d-wave superconductors 
[19]- _ 

Similarly, the ratio K c (T)/n a (T) can be computed and compared to the experiments. 
Within our model, it is given by 



K c (T)/ K a (T) = 0.2703 
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which is shown in Fig. 3 along with the thermal conductivity measurements of Ref. [11]. 

These expressions give a very reasonable description of the thermal conductivity for 
T/T c < 0.3. We note that a similarly good description of the thermal conductivity is given 
by the hybrid gap proposed in Ref. [11]. At higher temperatures, T/T c > 0.3, our simple 
model fails to describe the measured thermal conductivity, possibly due to the fact that 
phonons begin to play an important role as we approach T — > T c . Nevertheless, we can 
conclude that chiral d-wave SC is consistent with the experimental data of Refs. [11, 12] in 
the relevant low-temperature regime. Note also, that our calculations predict an interesting 
upturn in the ratio k c (T) / n a (T) as the temperature is further lowered. This prediction can 
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Fig. 3: Ratio of thermal conductivities of the <f||c and q\\a direction, plotted as a function of T/T c . 
The symbols represent the experimental data from Ref. [11]. 

be scrutinized experimentally, and may serve as a means to distinguish the present theory 
from the hybrid gap model that was proposed earlier. 

In the present context, the unconventional superconducting order in CcRhIn5 is of great 
interest. Let us briefly contemplate on the doped case. Inspecting Fig. 3 of Ohira-Kawamura 
et al [12] we may conclude that the order parameter in CeRhi_ x Co x In5 should be d-wave SC 
with an angular dependence / = cos(20), whereas the order parameter in C eRh\- x Ir x ln% is 
consistent with chiral d-wave superconductivity with an angular dependence / = e ±l< ^ cos(x). 
Therefore, the above approach will provide a basis to identify the many competing phases 
of the 115 compounds. Also, the phase diagrams for CeRh 1 _ x Co x In 5 and CeRhi_ x Ir x In5 in 
Ref. [12] are of great interest for the perspective of the Gossamer superconductivity, i.e. a 
phase with competing order parameters [10,20,21]. We observe that (a) the incommensurate 
phases in both CeRhi_ x Co x In5 and CcRhi_ x Ir x In5 are conventional spin-density waves, (b) 
the commensurate phase in CeRhi_ x Co x In5 and the incommensurate+commensurate phase 
in CeRh 1 _ x Ir x In 5 have d-wave symmetry. Therefore, there is a wide region where d-wave 
superconductivity coexists with unconventional nodal spin density wave order. 

In summary, we have successfully applied a nodal weak-coupling BCS theory to fit recent 
experimental data on the directional thermal conductivity of CeRhlns. We find that in 
contrast to CeCoIns, which has plain d-wave order, this compound is consistent with chiral 
d-wave superconductivity. Furthermore, this technique will allow us to identify the many 
different phases which were recently discovered in doped derivatives of these materials. 

* * * 

We would like to thank Y. Matsuda for very informative discussions on the supercon- 
ductivity in Celrln 5 . 
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